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Abstract 

We report low temperature, microwave transmission measurements on a new switchable and tunable class of nonlinear 
metamaterials. A wireless two dimensional array of Josephson junctions (JJ) is probed as a metamaterial where each plaquette 
in the array is considered as a met a-atom. In the presence of microwaves, this compact metamaterial of 30,000 connected 
meta-atoms synchronizes the flow of Cooper pairs to yield a single robust resonant signal with a quality factor of 2800 at the 
lowest temperature for our measurements. The transmission signal is switched on and off and its amplitude and frequency 
are tuned with either temperature, incident rf power or dc magnetic held. Surprisingly, increasing the incident rf power above 
a threshold causes the resonance to split into multiple discrete resonances that extend over a range of 240 MHz for a wide 
temperature window. We posit that this effect is a new incarnation of the inverse ac Josephson effect where the Josephson 
plasma frequency locks to the rf drive frequency generating quantized frequency steps, instead of voltage steps, in the absence 
of a dc bias. TgX . 


Prompted by theory im extensive studies of care¬ 
fully designed artificial materials, known as metamate¬ 
rials, are underway to find new electromagnetic phe¬ 
nomena such as amplification of evanescent waves |3j 
- [9], extraordinary transmission |10] , [II], negative 
index of refraction [12], [13], and classical analogues 
of electromagnetically-induced transparency m - eg 
However, in spite of recent advances using metallic meta¬ 
materials, requisites for practical applications such as 
switching on and off a signal and modulating a signal’s 
frequency and amplitude using external control param¬ 
eters still need to be addressed. Since fulfilling these 
conditions with existing metallic structures is challeng¬ 
ing, because of absorption losses, developing a new class 
of tunable metamaterials, such as superconducting meta¬ 
materials, is arguably warranted. m-m 

In this Letter, we show that tunable metamaterials 
already exist in the form of superconducting Josephson 
junction, JJ, arrays. Most remarkable, however, isn’t 
their tunability per se , but discovering an unanticipated 
rf power and temperature dependent manifestation of 
discrete frequency steps that, to our knowledge, has not 
been previously reported. This physical phenomenon 
opens up new possibilities for signal processing if scaled 
to THz frequencies m, new opportunities for high fre¬ 
quency detectors for radio telescopes [28] , and raises new 
theoretical questions about the relevance of the inverse ac 
Josephson effect to microwave transmission experiments. 
The inverse ac Josephson effect is the reverse of the ac 
Josephson effect. In the ac Josephson effect, an applied 
dc voltage across a JJ generates a periodically oscillat¬ 
ing supercurrent, pairs of electrons, that emits electro¬ 
magnetic radiation whereas in the inverse ac Josephson 
effect, shining electromagnetic radiation onto an unbi¬ 
ased JJ induces a dc voltage difference across the tunnel¬ 
ing barrier. Here we show that shining microwaves on a 
wireless Josephson junction array induces seven discrete 



FIG. 1: (Color online) Experimental set-up and sample a.) 
Schematic of the experimental setup containing a network an¬ 
alyzer and a waveguide, b.) Orientation of sample to electro¬ 
magnetic waves in waveguide, c.) SEM images of the JJ meta¬ 
atom: above, image of the tunnel junction which is shown in 
the circular indentation in the image and below, two images 
of a single JJ meta-atom with 4 JJ per loop, left image, in 
black and white, and right image, in color, d.) SEM image of 
the JJ array showing only a small fraction, ~80 meta-atoms, 
of the sample. 


frequency steps instead of voltage steps, and in between 
these resonances, there are narrow opaque windows in the 
transmission line. The inverse ac Josephson effect only 
applies for underdamped JJs [29, which also describes the 
JJs studied here. 

Two dimensional Josephson junction arrays are macro¬ 
scopic quantum materials with a highly tunable super¬ 
conducting current. We consider a wireless two dimen¬ 
sional network of 30,000 connected Josephson junction 
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(JJ) plaquettes with each plaquette having a square loop 
geometry containing 4 equally spaced JJs; we consider 
these plaquettes as meta-atoms which form, in aggregate, 
a switchable and tunable metamaterial. In contrast to re¬ 
cent microwave transmission reports on discrete split ring 
JJ meta-atoms m , eb , we find a sharp, robust reso¬ 
nance. As with our previous measurements on this sam¬ 
ple, abet at evanescent frequencies [9], the resonance is 
clearly distinct from noise without processing the data or 
using any additional microwave components, such as am¬ 
plifiers, attenuators, directional couplers or filters, which 
also differentiates our work from others. 

The experimental apparatus, described in detail else¬ 
where [9j, consists of a vector network analyzer connected 
to a X-band waveguide via phase-maintaining co-axial ca¬ 
bles as seen in the schematic of Figure 1 a. The waveg¬ 
uide, enclosed in a vacuum can, is cooled to cryogenic 
temperatures; it is magnetically shielded by (i— metal to 
reduce the ambient field. The JJ metamaterial is posi¬ 
tioned in the center of the waveguide and oriented such 
that the rf B-field is normal to the plane of the array as 
seen in the schematic in Figure 1 b. 

The sample is a portion of a much larger array [52] , 
[33]; its transport properties, before removing electri¬ 
cal contact pads and reducing its size, have been ex¬ 
tensively studied elsewhere. [34] Scanning electron mi¬ 
croscope, SEM, images of a single tunneling junction, 
Nb-amphorous silicon-Nb, and of a single meta-atom are 
seen in the upper and lower left images, respectively, of 
Fig. 1 c. The 500 nm niobium thick square loop is 29 
/im x 29 iim with a linewidth of 5 nm. Each of the 4 JJs 
in the loop has an area of 2.5 /im 2 with an insulating 
barrier of 10 nm. A colored optical microscope image of 
a single meta-atom showing niobium coverage is seen in 
the lower right image of Fig. 1 c; the non-continous cov¬ 
erage of niobium indicates that resonance occurs because 
charged particles must tunnel across capacitive tunneling 
barriers to circulate throughout the network. A portion 
of the JJ metamaterial is seen in the SEM image of Fig. 
1 d. The individual JJ meta-atoms are about 1500 times 
smaller than the microwave wavelength, which is approxi¬ 
mately 40 mm, placing them in the deep sub-wavelength 
limit. j35l The superconducting transition temperature, 
T c , of the individual niobium islands is 9.2 K. [36! 

If we consider that each JJ contributes a capacitance, 
C, estimated to be 0.9 pF, a normal state resistance, R, 
estimated to be 510 12, and a critical current, 7 C , which 
is estimated to be 0.165 /iA, m we can estimate the 
degree of damping in a junction, using the dimensionless 
McCumber parameter, /3 C , which is (3 C = 2eIc ^ c ~ 118, 
where e is the electron charge and h is Planck’s constant 
divided by 27r. Since (3 C > 1, the junction is said to be 
underdamped meaning that there are inertial effects of 
the ’’phase” particle [38] which satisfy one of the condi¬ 
tions required for the inverse ac Josephson effect [29]. A 
second condition that needs to be met is that the driv- 
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FIG. 2: (Color online) Switching a resonance on and off with 
temperature and rf power, a.) Temperature dependence at 
fixed rf power, -55 dBm, for T = 2 K (red) and 10 K (black) 
and b.) rf power dependence at fixed temperature, 2 K, for 
rf power -55 dBm (red) and -5 dBm (black). Note that the 
range in transmission on the vertical axis is different in a 
and b. Also note that the transmission values above 0 dB 
correspond to 0 dB: this offset is due to the calibration being 
done at room temperature while the measurements are at 
liquid helium temperatures. 


ing rf frequency should be significantly larger than the 
natural Josephson plasma frequency, / = ~ 15 

GHz. [39] The resonant frequency, which is the plasma 
frequency [9], also depends, in part, on the sample’s loca¬ 
tion in the waveguide indicating the need for additional 
capacitive terms [40 , which would lower the resonant 
frequency but not by enough to satisfy this second con¬ 
dition. Nevertheless, we will show that the mechanism 
responsible for the transmission features seen here are 
due to the inverse ac Josephson effect. 

Cooling the sample to temperatures below the su- 
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FIG. 3: (Color online) Tuning the resonance with temper¬ 
ature and rf power. Temperature dependence when the rf 
power is fixed at a. )-55 dBm and b.) -49 dBm for tempera¬ 
tures: T = 2.0 K(black), T = 3.0 K (red), T 4.0 K(green), T 
= 5.0 K (blue)= 6.0 K (purple), and T = 7.0 K (brown). Ag¬ 
gregation of temperature-dependent data at fixed power: c.) 
-55 dBm and d.) - 49 dBm. Each color represents a different 
temperature. 


perconducting transition temperature and applying low 
power microwaves to the JJ array results in a sharp reso¬ 
nance at f — 8.08 GHz appearing in transmission as seen 
in the plots of Fig. 2. Extinguishing the resonance is 
achieved by either increasing the temperature above the 
superconducting transition temperature as seen in Fig. 2 
a, or by increasing input power as seen in Fig. 2 b. 

Increasing temperature, at fixed rf power, shifts the 
resonance to lower frequencies and broadens the reso¬ 
nance. This is seen in the transmission vs frequency plots 
for different temperatures in Fig 3 a-d. As temperature 
increases, Cooper pairs break lowering the critical cur¬ 
rent, I c , and thus shifting the resonance to lower frequen¬ 
cies. While the density of Cooper pairs decreases, the 
density of normal conduction electrons increases, lead¬ 
ing to more dissipation and insertion loss; for instance, 
a jump of - 5 dB in insertion loss is measured when the 
temperature reaches 7 K as seen at 7.2 GHz in Fig. 3a 
and 3b. While these behaviors are reminiscent of those 
seen with superconducting metamaterials [22 , [42] , they 
are not identical. Here, the resonances are more sensitive 
to temperature increases due to Cooper pairs tunneling 
across barriers as described by the Ambegaokar-Baratoff 
current-temperature relation. |43j 

Similar temperature dependent resonant responses 
persist when comparing responses at two different rf pow¬ 
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FIG. 4: (Color online) RF power dependence at fixed temper¬ 
ature, T=l.l K. The black arrows indicate multiple discrete 
resonances that appear below the central resonance and at 
high rf power; these are spaced about 0.04 GHz apart. Inset: 
RF power dependence at fixed temperature, T = 1.24 K for 
-49.5 dBm (blue) and -49 dBm (red). 


ers, -55 dBm, as seen in Fig. 3 a, and -49 dBm, as seen 
in Fig. 3b. At -49 dBm, the higher rf power, the res¬ 
onance also shifts to lower frequencies while broadening 
with increasing temperature; however, the overall ampli¬ 
tude decreases and the resonance dip splits into two for 
T < 6 K. 

To elucidate the effect of the resonances splitting at - 
49 dBm, we compare two large sets of temperature data, 
one at -55 dBm and the other at -49 dBm as seen in 
Fig. 3 c and 3 d. The aggregate temperature data at 
-49 dBm and for T < 5.2 K clearly show a simple ”egg 
crate” pattern in transmission that is not present at -55 
dBm. With increasing rf power, screening currents begin 
to slosh back and forth generating an ac magnetic field 
which cancels the applied RF field. Furthermore, when 
these induced ac supercurrents, driven by external mi¬ 
crowave radiation, lock to the natural oscillations of the 
array there is a range of DC currents for which there are 
constant frequency steps. This current range is depicted 
in the number of resonances, with different amplitudes, 
appearing at the same frequency. 

Frequency locking is also manifested when increasing 
the RF power at fixed temperature, T= 1.1 K, as seen 
in the transmission vs frequency plot in Fig. 4. As the 
RF power increases, the number of dips also increase. 
While an increase in the current decreases the depth of 
the dip, it does not explain multiple dips. We surmise 
that the increase in the number of dips is analogous to 
Shapiro steps observed in dc current-biased experiments 
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FIG. 5: (Color online) a.) Resonant frequency versus temper¬ 
ature at - 49 dBm. Seven steps are clearly visible for T < 5 K. 
Each of the different colors correspond to separate dips con¬ 
taining multiple resonances. A line is draw to show the linear 
dependence at temperatures close to the transition tempera¬ 
ture of the array. Inset: Comparison between -49 dBm (red) 
and - 55 dBm (blue), b.) Transmission Amplitude, S 21 , ver¬ 
sus temperature at -49 dBm (red and green) and - 55 dBm 
(blue). 


with a microwave irradiated JJ. [431 Here, however, there 
is no applied bias, only an induced dc voltage bias. Since 
the device is nonlinear, locking occurs not only at ojrf, 
but also at harmonics of the RF frequency - that is, 
at multiples of ujrf , such that twjrf = 00 j - giving rise 
to discrete dips appearing at lower frequencies and sep¬ 
arated by 0.04 GHz as indicated by the black arrows in 
Fig. 4. A close-up image of a resonance splitting into 
two with higher RF power for T =1.24 K is seen in the 
inset of Fig. 4. 

We also analyze these results by studying how the reso¬ 
nant frequency and transmission amplitude change with 
temperature and rf power. These results are shown in 
the plots of Fig. 5. There are seven steps in resonant 
frequency below 5 K that are clearly observed, with each 


step represented by a different color, in Fig. 5 a. Con¬ 
sidering the temperature to be proportional to the dc 
current and the frequency proportional to dc voltage, by 
the Josephson relation, Vr>c = n{hjujRF)/ 2e , then these 
steps look similar to those found in I(V) curves for spon¬ 
taneous emission in biased, underdamped JJ arrays. [31] 
The range of phase locking and DC currents in our experi¬ 
ments is indicated by the plateaus for each frequency; the 
width of the plateau becomes smaller with higher order 
steps because the current is expressed in terms of a Bessel 
function. [38] Between resonant frequencies, the phase is 
unlocked leaving narrow electromagnetically opaque re¬ 
gions. A diagonal line is draw in the plot to show the 
linear behavior between frequency and temperature at 
temperatures close to T co , the transition temperature of 
the array,~6.5 K. Comparing resonant frequencies for -55 
dBm and -49 dBm reveals two main differences: first is 
the appearance of steps at -49 dBm and second, is the 
systematic higher shift in the frequency vs temperature 
at -55 dBm below the transition temperature as shown 
in the lower inset of Fig. 5 a. This shift is due to the 
induced rf currents being lower. 

The strength of the resonance, transmission amplitude 
of S 21 , is dependent on rf power as seen in Fig. 5 b. 
However, we find that it also varies with temperature, in 
contrast to split ring resonantors with JJs m, except 
at temperatures between ~5 K and 6 K which coincides 
to temperatures where frequency steps no longer prevail. 
Since the resonance splits at -49 dBm, we label the two 
resonances either left as seen in green in Fig. 5 b or 
right as seen in red in Fig. 5 b. The -55 dBm and - 
49 dBm curves crossover at ~6.5 K which is the same 
temperature that the resonant frequency flattens out as 
seen in Fig. 5 a. This is indicated by a vertical dash line 
in Fig. 5. There is still a resonance beyond 6.5 K that 
doesn’t change in frequency suggesting superconducting 
quasiparticles tunneling until T ~ 8.5 K. 

In conclusion, our results demonstrate that a wire¬ 
less JJ array in a waveguide can controllably switch on 
and off a microwave signal with temperature and input 
power. The array can also modulate the signal over a 
range of frequencies and amplitudes by varying either in¬ 
put power, temperature, or dc magnetic field (not shown 
here). Modulating signals with temperature is tradition¬ 
ally considered to be quite slow; however, for JJ arrays, 
the response occurs rapidly and without noticeable de¬ 
lay. The signal is sharp and robust without processing 
the data and does not require the aid of amplifiers and 
filters in the experimental configuration. Moreover, our 
results reveal qualitatively new phenomena for tunable 
metamaterials that occur as the input power increases, 
breaking the central resonance down to discrete frequen¬ 
cies. Finally, this alternative perspective on metamate¬ 
rials may provide new strategies for future metamaterial 
design and tunability. 
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